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We show how a single, harmonically trapped atom in a tailored magnetic field can be used for 
simulating the effects of a broad class of non-abelian gauge potentials. We demonstrate how to 
implement Rashba or Linear-Dresselhaus couplings, or observe Zitterbewegung of a Dirac particle. 



Simulating complex quantum systems with the help 
of trapped cold atoms has become a flourishing branch 
of quantum optics. The exquisite control and complex 
internal structure of cold atoms have allowed successful 
experimental simulation of systems and effects ranging 
from quantum phase transitions in the Bose-Hubbard 
model over Anderson localization 0, Q, to cosmo- 
logical models [3]. Recently, a lot of research efforts 
have been directed towards realization of Abelian or Non- 
Abelian gauge fields. By bathing an optical lattice in ad- 
ditional weak, non-resonant light that can create Raman- 
transitions between hyperfine levels, one can create ar- 
tificial magnetic fields which can be extremely strong. 
Vortex formation in a BEC has been observed due to 
such artificial magnetic fields Even the creation of 
magnetic monopoles has been proposed @. Non-abelian 
gauge fields can be useful for studying spin-tronics mate- 
rials with various spin-orbit couplings, Berry phases, or 
topologically protected qubits. 

Most of these quantum simulations have been proposed 
or performed for atoms trapped in an optical lattice or 
Bose-Einstein condensates, requiring rather sophisticated 
experimental setups. Many physical systems are, how- 
ever, interesting as single-particle systems. Such is the 
case e.g. for Anderson localization, the low-energy be- 
havior of electrons in Graphene, or the relativistic mo- 
tion of electrons that leads to the effect of Zitterbewe- 
gung 0. It would be highly desirable to have a simple, 
versatile system that allows to simulate such single par- 
ticle dynamics. We show here that such a system can be 
constructed from a single, harmonically trapped atom ex- 
posed to a suitably tailored real (physical) magnetic field. 
We show that by simply changing some gradients of dif- 
ferent field components different non- Abelian gauge fields 
can be simulated, giving rise, for example, to Rashba or 
linear Dresselhaus coupling. We also show that Zitterbe- 
wegung should be easily observable. Abelian gauge fields 
can also be obtained, but are less interesting in the pro- 
posed setup, as they do not depend on position. 

Consider a single atom of mass m trapped in a har- 
monic potential (frequency uj) and exposed to a magnetic 
field B(x, t). For concreteness, let us assume a neutral 
atom in an optical dipole trap. The magnetic field should 
not be used for any part of the trapping, but be control- 
lable independently of the trapping. We will also assume 
that the atom is cooled close to the ground state such 



that approximating the trap by a harmonic potential in 
all three directions is indeed reasonable. The hamiltonian 
of the system reads 
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gJ^iB'B{^,t)■3/h, (1) 



where gj, ^b, and J are the g-factor, Bohr-magneton, 
and total angular momentum of the atom, and we have 
neglected the nuclear spin contribution in the Zeeman 
term. Suppose that the magnetic field varies slowly on 
the scale of the trapping potential, such that we can ex- 
pand it in a power series about the origin x = 0, 



B(x,i) =B(0,t)+^VB., •xej + e'(x2), 
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with unit vectors in directions i = 1,2,3 = x,y,z. 
We have neglected higher order terms starting with the 
second order. One can still include the second order, with 
consequences to be discussed below, but for the moment 
suppose that the magnetic field varies slowly enough over 
the length-scale of the atomic motion in the trap that 
the quadratic term in x from expanding B can indeed be 
neglected compared to the quadratic term describing the 
trapping potential. 

The linear term in x then leads to a coupling x • V (B • 
J) = x • J^iC^^i lx=o)"^i- ^ow let us canonically trans- 
form x/xq — >■ — p/poi p/po ~^ x/iq, where I have in- 
troduced the natural scales xq and po of the canoni- 
cal coordinates and momenta of the atom in the trap, 
xq = y^h/{muj) and po = y/hmtu. The harmonic os- 
cillator part in the hamiltonian is invariant under this 
transformation, but the x in the Zeeman term becomes 
a — p, and we get the transformed hamiltonian H, 
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H^-^ + -mcj V - -p . A + gJ^lBB{0, t) ■ 3/h , (3) 
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Note that A is independent of position and thus com- 
mutes with p. We can therefore rewrite the hamiltonian 
in the form 

i/^ -^(p-eA)2 + imc^2^2^0+^B5jB(O,t)-J/ft, (5) 
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which makes appear A as a gauge potential. The con- 
stant (j) s-cts only in the angular-momentum Hilbert 
space, 
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and can be considered an anisotropy for the angular mo- 
mentum. We can write it as 
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where the inverse tensor of inertia / ^ has matrix ele- 
ments 
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The effect created by the gauge potential plays an appre- 
ciable role only if the components of eA arc comparable 
to the corresponding components of p. The latter are, 
close to the ground-state of the harmonic potential, of 
the order of po. With J ~ S, we are thus led to a con- 
dition for the magnitude of the magnetic field gradient 
B' 
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which should be satisfied for at least one component 
of B. Inserting the data for ®''Rb, we find a gradient 
B' — 2.1Gauss/mm for uj = 27r-lkHz, which appears to 
be a very convenient order of magnitude. 

It is worthwhile spelling out the components of A ex- 
plicitly. Up to the common prefactor ^^^f^- we have 

A^ cx dxB^ + dxBy Jy + dxB^ (10) 

Ay cx dyBx Jx + dyBy Jy + dyB^ (ll) 

A, cx dAJx + d,ByJy + d,B,J,. (12) 

The components of J are real physical angular momen- 
tum, and thus satisfy [Jx^ Jy] = ihJz in any Hilbert space. 
Therefore, obviously, the gauge field is in general non- 
abelian. The remarkable thing about A is that it is eas- 
ily programmable by an appropriate choice of magnetic 
field gradients. Choosing for instance only a field gradi- 
ent of one magnetic field component leaves us with only 
one component of the angular momentum and we thus 
get an abelian gauge potential. Choosing arbitrary lin- 
ear combinations of angular momentum operators is only 
restricted by the properties of the real physical magnetic 
field employed, in particular divB = 0, as we shall see 
below. 

Before working out a few specific examples, three more 
remarks: 

1. We chose the same trap frequency w in all direc- 
tions, which is certainly somewhat unusual for an 



optical dipole trap. After canonical transformation, 
this translates into the same mass for all three di- 
rections. Had we different trap frequencies in dif- 
ferent directions, we would get effectively different 
masses in different directions, i.e. a (diagonal) ef- 
fective mass tensor, similar to the situation in a 
semiconductor. This may or may not be useful, 
depending on what one wants to simulate. 

2. Similarly, if we kept the quadratic term in the ex- 
pansion of B, we would get additional quadratic 
terms (which would also depend on the internal 
state of the atom) contributing to the trapping po- 
tential, such that, again, after the canonical trans- 
formation we would end up with different masses 
in different directions. 

3. What was derived here for a single atom should 
translate immediately to a BEC, if the interactions 
between atoms are negligible. If they are not, any 
interaction potential y(xj— x^) will translate into a 
momentum dependent interaction after the canon- 
ical transformation. 

Let us now look at some concrete examples of non- 
Abelian gauge potentials, inspired by the list in I 
will assume that we have an atom with J 1/2, so that 
the Ji are given by Pauli- matrices, Ji = {h/2)ai. 

Rashba coupling. This is a coupling a{py<7z — Pzay), 
with some constant a. It can easily be achieved by choos- 
ing dxBi = 0, for i = x,y,z, dyBy = = dzB^, and 
dyBz = B' ^ dzBy. No magnetic mono-pole is required, 
divB = 0, but the field has a finite curl in x-direction, 
(V X 'B)x = 2B'. That curl is needed at x = 0, where 
we did the expansion of the magnetic field. But since we 
insisted on a linearly growing field over the lengthscale 
of the motion in the trap, we request basically circular 
magnetic field lines around the x-axis, with magnetic field 
strength growing linearly with distance from the x-axis. 
That is, of course, quite different from what a current car- 
rying wire along the x— axis would do (apart from the fact 
that that wire would have to go straight through x = 0, 
where we would like to trap the atom) . A possible way of 
creating the required magnetic field might be through an 
electric field that increases linearly in time by using seg- 
mented cylindrically arranged electrodes. Choosing the 
right voltage profile for all the electrodes one can create 
a radially symmetric electric field oriented in x-direction 
that increases proportionally in radial direction, and gen- 
erates the desired B-field according to Maxwell's equa- 
tion V X H = i) -I- j. 

The inverse tensor of inertia is here simply 
I^^ = diag(0, i?'^, B'^) and just leads to a constant 
(j) — B'^1, i.e. an identity matrix in spin Hilbert-space. 



3 



Linear Dresselhaus coupling. Here we want Ay = aay, 
Az = — cttTz- Choose again dxBi = Vi, but then 
dyBy = _B' = —dzBz, and all other derivatives equal 
zero. Again, this is compatible with divB = 0, and we 
basically get a field that is point symmetric, growing lin- 
early in radial direction, but pointing towards the center 
on the z— axis and outward on the ?/— axis, like from a 
quadrupole magnet. 

Inverse tensor of inertia and the resulting potential (p 
are identical to the Rashba case. 

Graphene sheet in vicinity of Dirac point. This case is 
more problematic: we want Ay = aay, Az = aaz, which 
leads to dyBy = B' = dzBz. But divB = then requires 
dxBx = —2B', and thus adds an extra coupling term 
~2pxOx to the desired Py(Jy -\- PzCz- 

Zitterbewegung. Zitterbewegung (ZB) is an interference 
effect first predicted by Schrodinger for relativistic 
spin- 1/2 particles that leads to a jittering motion on the 
length scale of the Compton wavelength of the particle, 
i.e. h/mc ~ 10~^^m for an electron 0. This short 
length scale has so far prevented direct experimental 
observation of ZB for relativistic electrons. However, 
the effect should exist for any spinor system with linear 
dispersion relation. Consequently, ZB has been studied 
theoretically in several systems, including mesoscopic 
wires 0, graphene [l^l, ion traps 11 1, and optical 
lattices [8|. Rashba couplings in quantum dots with 
rotationally invariant potentials in 2D was studied in 
[l^ . Very recently, ZB was observed experimentally with 
a single trapped ion [l3| . ZB has mostly been studied for 
a free particle. The ZB then dies out after a short time, 
when the two wave-packages corrcsonding to spin-up 
and spin-down have separated enough to prevent further 
interference. It is interesting to consider how the addi- 
tional harmonic confinement potential in ([5]) modifies 
the ZB. One might expect that the confinement will in- 
crease the time interval in which the ZB can be observed. 

We now study ZB for the harmonic confinement ([5|) 
in the Rashba case. We show that ZB exists even if the 
harmonic oscillator is initially uncxcited, and that it can 
persist for arbitrarily long times. Note that the third 
direction, x in the above notation, is not affected by the 
gauge field and just separates. It is convenient then, to 
rewrite ([5|) in terms of annihilation (creation) operators 
a(a^) and b{h^) for the y and z components, respectively. 
We are thus lead to 

H = hw[a}a + b'^b + ibi[- {a- a))az + {b -b'^)ay) 
+bo • 0-) , (13) 

where the dinicnsionlcss parameters arc 61 = Q-jf^sB 

and bo = , and we have suppressed an irrelevant 

constant. 



It is straightforward to express H in basis states \nms) , 
where n^m € {0, 1,2,...} are the occupation numbers for 
oscillators a and b, and s £ {±1} label eigenstates of cr^ 
corresponding to these eigenvalues, and to diagonalize 
the Hamiltonian numerically. This basis turns out to be 
highly suitable — taking into account only up to 5 exci- 
tations per oscillator (i.e. a 72 dimensional basis) allows 
one to find the lowest 26 eigenstates already with 4 sig- 
nificant digits for bi =0.1 (compared to the case with 10 
excitations per oscillator, or 242 basis states, which we 
used in the numerical simulations). With the obtained 
propagator wc can study the time evolution of the aver- 
ages of the observables {y{t)), {py{t)) etc. 

Figure [T] shows {y{t)) for bo = 0, different values 
of bi, and the two initial states \ip{0)) = |00+) and 
|00— ). Both {z{t)) and {pz{t)) remain always zero in 
these cases, whereas (y(i)) shows oscillations. They 
are always with opposite sign for these two different 
initial states. For small values of bi (5i = 0.1) the 
oscillations appear periodic, whereas with increasing 
61, additional harmonics appear that make the signal 
look more and more erratic. By changing the initial 
state of the spin, the direction of the oscillation can 
be chosen. E.g. for|7/)(0)) = (|00+) +i\00-))/V2, we 
obtain {y{t)) = = {py{t)), whereas the z-componcnt 
now shows the signal we had for {y{t)). The state 
IV'(O)) = (cos(7r/8)|00-H) + isin(7r/8)|00-)) leads to 
oscillations with {y{t)) = —{z{t)). It appears thus that 
the ZBis always one-dimensional as long as the initial 
state is chosen in the subspace {|00-f), |00— )}. The 
superposition |V'(0)) = (|00-h) + |00-))/V^ switches the 
ZB off in both components. Note that one-dimensional 
ZB was also predicted for a ID harmonic confinement 
and a specific initial spin state, with the ZB perpendic- 
ular to the free ID motion Q. 

The results for small bi can be easily understood ana- 
lytically, by going to the Heisenberg picture and expand- 
ing the time dependent operators to lowest order in bi. 
We find, correct to order 0(61), at bo = 

{y{t)) ~ (y(0))cosa;t+ (py(0))sina;t 

+\/2fei ((72(0)) sin (14) 
{z{t)) ~ (z(0)) coswi + (^^(0)) sinwt 

-%/2foi(cry(0)) sinwt (15) 
{py{t)) ~ {py{0)) COS ujt — {y{0)) sinujt 

+^bi{az{0)){cosut-\) (16) 
{pz{t)) (z(0))coswi- (z(0))sintjt 

-V2bi{ay{Q)){cosLot-l) . (17) 

All positions and momenta are expressed in terms of the 
length-scales xq and momentum scales po of the harmonic 
oscillator, respectively. The corrections to Eqs. (|14lir7|) 
are of order 0{B\). We recognize in the terms indepen- 
dent of 61 the unperturbed motion of the 2D harmonic 
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FIG. 1: (Color online) Zitterhewegung in {y{t)) (in units of 
harmonic oscillator length xq) for initial states |00+) (full 
lines) and |00— ) (dashed lines); bo = 0, and bi takes the 
values bi — 0.1, 0.5 and 1 for black lines, green (light grey) 
lines, and blue lines (dark grey), respectively. 



FIG. 2: (Color online) Trajectory of 2D Zitterhewegung in 
{{y{t)) , {z{t))) (in units of harmonic oscillator length xq) for 
an initial state |a/3+), where |a) and \(3) are coherent states 
with a = 1 and /3 = i at bo = and bi = 1. The trajectory 
is shown for < ojt < 50. The initial state is localized at 
(j/(0)) = x/2 = (p40)), (z(0)) = = {py(0)) and would lead 
for fei = to circular motion in the yz-plane. 



oscillator. This motion can be switched off by chosing 
average initial positions and momenta equal zero, as is 
the case for the initial states discussed above with the 
two harmonic oscillators initially in their ground state. 
All motion is then entirely due to the ZB and indeed ID, 
in a direction given by the vector ((^^(O)), ((Ty(0))). In- 
terestingly, at |&i| ^ 1 the ZB is itself harmonic with 
an amplitude controlled by 6i. The fact that the length 
and momenta scales of the ZB are set by the harmonic 
confinement should allow for a much simpler experimen- 
tal verification than for Dirac electrons. Moreover, the 
harmonic confinement potential keeps the wavepackages 
corresponding to the two different spin components to- 
gether, preventing the decay of the ZB, which should fa- 
cilitate its experimental study. A damping of the ZB can 
arise at higher values of bi and initial coherent states of 
the harmonic oscillators, which will in general be smeared 
out due to the anharmonicity mediated by the coupling 
to the spin (see Figsl^] and At & = 1 the dynam- 
ics looks random and diffusive in the j/z-plane, but the 
trajectory of average values {{y{t)), {z{t))) is of course 
entirely deterministic and reproducible. 

The increasingly erratic behavior of the ZB with in- 
creasing bi can be understood by looking at the spec- 
trum of H, plotted in Fig|4l The oscillator states at 
&i = containing n energy quanta which are 2{n -\- 1)- 
fold degenerate (states |nO-(-), |n — 1 1+), . . . , |0n+), and 
the same set once more but with spin down). They split 
with increasing bi and lead to several avoided crossings. 
The resulting incommensurate frequencies lead to the ob- 
served quasi-periodic, apparently random behavior. 

In summary, we have shown how a single trapped atom 
with hyperfine structure trapped in a magnetic field with 
suitably tailored field gradients can be used to simu- 
late the effect of non-abelian gauge potentials. We have 
demonstrated how different effective spin-orbit couplings 
(such as Rashba or linear Dresselhaus couplings) can be 
easily obtained, and we have proposed a new way of ob- 




FIG. 3: (Color online) Probability distribution p{y, z) of find- 
ing the atom localized at position y,z (in units of xo) at 
ujt — 10, for bi = 1.0, b = 0. The initial coherent state (same 
state as in Fig[2]) gets smeared out due to the anharmonicity 
created by the Rashba coupling. 



serving Zitterhewegung of a harmonically trapped parti- 
cle. An immediate consequence of the fact that ([S]) was 
obtained by a canonical transformation that exchanges 
position and momentum, is that the ZB in y(t) will, of 
course, show up in Py{t) in the original system. The 
simplicity and flexibility of the proposed setup may also 
allow the study of other spin-orbit couplings, as well as 
of applications such as robust quantum gates based on 
topological phases [14 \ . 
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